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THE DEVIATIONS OF FALLING BODIES. 

Bt F. R. Motjlton.* 

1. Introduction. The subject of the deviations of freely falling bodies 
has been considered in many memoirs from the time of Laplace and Gauss 
to the recent work of Roever and Woodward, f All writers have agreed 
that a body falling from rest near the surface of the earth wiU deviate to 
the eastward with respect to a plimib-line hung from the initial point, but 
a great variety of results have been obtained regarding the deviation 
measured along the meridian. For example, Laplace found no meridional 
deviation. Gauss a small deviation toward the equator, Roever a deviation 
toward the equator several times that of Gauss, and Woodward a small 
deviation away from the equator. The diversity of these results is the 
occasion for this paper. 

In the present paper the solution of the problem of the deviations of 
falling bodies is made to depend upon two well-estabhshed mathematical 
theories, viz., that of the solution of analytic differential equations in the 
neighborhood of regular points, and that of implicit analytic fimctions. 
The methods employed are rigorous, of general appUcabiUty, and relatively 
simple. The problem involves three parameters, viz., the height from which 
the body falls, the rate of rotation of the rotating body, and the oblateness 
of the rotating body. Only the first of these is arbitrary, and it is shown 
that the deviations of the falling body can be developed as converging power 
series in it if it is sufl&ciently small. In the meridional deviation the 
coefficient of the first power of this parameter, which is a function of the 
other two parameters and can be developed as a power series in them, is 
identically zero, at least in the case where the rotating body is a figure of 
revolution; but the coefficient of the second power of the parameter is not 
zero. The sign of this coefficient rigorously determines the character of 
the deviation when the distance through which the body falls is not too 
great. The principal part of the coefficient of the second power of the 
parameter agrees exactly with the result found by Roever in his first paper 
by an altogether different method. Differences in terms of higher order 
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would be found if they were computed because his definitions of the devi- 
ations differ somewhat from those used in this paper. 

The question naturally arises why so great a variety of results has been 
obtained. The answer seems to he in the fact that the solution of the prob- 
lem depends upon two distinct mathematical processes, however much 
they may be concealed by the details of treatment; that both of them 
nattirally lead to power series in all three parameters; that the coefficients 
of all terms of the first degree in the first parameter are identically zero; 
and that methods of approximation are not adapted to establishing identi- 
ties. It is clear from what has been said that the correct result could 
never be established by expUcitly developing the results in aU three param- 
eters because only a finite number of terms could be computed; and it is 
equally clear that great care would be necessary to secure exactly the same 
approximation in solving the differential equations and the implicit fimc- 
tions by approximating methods because the two processes are so different. 
Most writers have developed the results from the beginning in all three 
parameters, and have used approximations without proving their legitimacy. 

Z» The Differential Equations. Take the origin at the center of the 
rotating body and the principal polar axis coincident with the axis of rota- 
tion. Then, if r represents the distance, (p the latitude, 8 = do + at + X 
(where w is the rate of rotation of the body) the longitude, and V the 
potential function, the differential equations of motion are 

(1) dtV^) + '[''+-dt) sm^cos^ = ^, 

The potential function F in all cases can be developed as an infinite 
series of spherical harmonics whose coefficients depend upon the distribution 
of mass of the rotating body. If the rotating body is a figure of revolution 
about the axis of rotation whose density does not depend upon the longitude, 
the function V can be developed as a series of zonal harmonics of the form 

(2) 7=^ + ^(l_3sin2^) + .... 

Only the first two terms have sensible effects upon f aUing bodies and plumb- 
lines for short distances. 

Let Po be the point from which the body falls and from which the 
plumb-line is suspended, Po'-"^ the point where the line from the origin 
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to Po pierces the surface of the rotating body, Pi the point where the 
falling body strikes the surface of the rotating body, and Pi the point where 
the plumb-bob touches the surface of the rotating body. For notation 

let the coordinates of Po be ro(l + K), <po, 0, 

let the co6rdinates of Po^°^ be ro, <po, 0, 

(3) 

let the coordinates of Pi be n, <pi, Xi, 

let the coordinates of Pt be r^, (pi, X2. 

The differential equations (1) are to be solved subject to the initial con- 
ditions 

r(0) = ro(l + h), ^(0) = <p„ X(0) = 0, 

r'(0) = 0, ^'(0) = 0, X'(0) = 0. 

The problem is to find the points Pi and Pi. If F is a sum of zonal 
harmonics, X2 = 0. When Pi and Pa have been found the longitudinal 
and meridional deviations in angular measure are respectively Xi — X2 
and <pi — <pi. If the first is positive the deviation is to the east, and if the 
second is positive the deviation is from the equator. 

3. Integration of Equations (1). Since equations (1) are regular in the 
vicinity of the initial conditions, it follows from the general theory of 
analytic differential equations that they can be integrated uniquely as 
power series in t which converge for \t\ sufficiently small. It follows 
from (4) that the solution has no terms of the first degree in t. Hence 
the solution in the general case in which F is a sum of spherical harmonics 
has the form 

r = ro(l + %) +0t + a^f + agf + ad" + • • •, 

(5) ^ = ^0 + Oi + htf + bs** + 64*^ + • • •, 

X = -I- Of + C2<2 + CzP + C4<* + • • •, 

where the Oj, hi, and c» are constants which can be obtained by substituting 

(5) in (1) and equating coefficients of corresponding powers of t. 

But if F is a sum of zonal harmonics, and consequently independent of X, 
it follows from the third equation of (1) that C2 = 0; and then from the 
first two equations of (1) that 03 = 63 = 0. Consequently, in this case the 
solution has the form 

r = ro(l +h) + aj? + Oi« + aj," + atV + • • •, 

(6) ^ = -Po + hit-' + 0<« + 64** + h,t' + • • •, 
X = + + csi' + C4f* + csf^ + • • •• 
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In general, the coefficients of all powers of t beyond the third are distinct 
from zero. 

The fiill details will be given only in the case where F is a sum of zonal 
harmonics. It is found that 



2a2 = wVo(l + h) cos Vo + {'df)> 



262 = — |w^ sin 2^0 — 



dV\ 



(7) 



6C3 = - 4C0 [,7(1^ - fc2 tan ^„ J , 
12a4 = 4ro(l + Kjbi' + oP cosVoOz - ro(l + ^)w^ sin 2<poh2 

+ 6ro(l + h) cosVoCa + \-q^^)^^^ + (.a^a^ j,^^' 

, „, 120262 <^ sin 2^0 o • o 2 «^„ o », 

12&4 = - ;:;(j^^ - If^^:^^ «2 - 3« sm 2^„cs - «^ cos 2^062 

The coefficients of higher powers of i are not needed. 

4. Determination of the Coordinates of Pi. Let i\ be the time at which 
the faUing body reaches the surface of the rotating body. Then the 
coordinates of Pi are obtained by replacing f by ii in equations (6). Hence 
the value of i\ must be found. 

The point Pi is subject to the condition that it shall be on the surface 
of the rotating body, which is an equipotential surface for the potential 
function Y and the rotational potential ^wV^ cos^ <p. Hence the coordinates 
of Pi satisfy the equation 

|wW cosVi +^ + |3(l - 3 sinVi) + • • • = C, 

where C is determined by any point on the surface. But the coordinates 
of Po^°^ also lie on the surface of the rotating body and satisfy an equation 
of the same form. The difference of these equations, which eliminates the 
CQustant C, is 

F = Win^aoB^ ^1 - u"- cos^ -po) + « (^ - ^) + /3 (^ - ;:^) 

(8) / • 2" • 2 \ 

If ii is put in (6) in place of t, and if the series for ri and (p\ are sub- 
stituted in (8), the resulting expression can be developed as a power series 
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in ti and h which will converge provided 1 <i | is sufficiently small and 
\h\ <1. That is, (8) takes the form 

(9) F = pQi, h) = 0, 

where p is a power series in h and ti. Moreover, it is easily foimd that 

(10) p(0, 0) = 0, ( 1^ )„ = 0, ( 1^ )„ = 4[(a,«'>)^ + ro^(b2«>)^], 

where 02*^ and 62'"^ are obtained from 02 and 62 by putting h equal to zero. 
It follows from the theory of implicit functions that (9) is solvable for h as 
a power series in V^, vanishing ior h = and converging for \h\ suffici- 
ently small, provided 4[(a2^"0^ + »'o^(&2^°0^] is distinct from zero. But this 
quantity is the square of the acceleration of gravity at the poiat Po^'\ 
and consequently is not zero. Therefore the solution of (9) has the form 

(11) h = V^pi(Va), 

where pi( V^) is a power series in V^. 

It has been stated that it is sufficient to have the developments of 
series (6) to terms of the fourth degree inclusive in t. If the series are 
terminated at this point and F is a sima of zonal harmonics, F becomes 
a function of (1'' and h, and equation (8) can be solved for ti^ as a power 
series in /i of the form 
(12) ti^ = aih + a2K'+ •••. 

The first two terms of this series, which alone are required, are identical 
with the first two terms of the square of (11). 

In order to find the coefficients ai and a2 the explicit form of (9) is 
reqidred. It is foimd from (6), (7), and (8) that 

F = Oioti^ + aoih + 020^1* + anti% + Ooih^ + • • • = 0, 
oio = 2[(o2«»)^ + r.'ih^o^n 
Ofli = 2ro02^°^ 
020 = 2c2(«a4<« + 2ro262^»^b4«> 



(13) 



+ [ W cos^ ^0 + ^ + ^5 (1 - 3 sm^ <po) ] (a2«'>)^ 

an = [ro«^ cosVo +^ + ^d - 3 sin^o) + 2 (^)J««^'^ 

+ [ - Vo)^ sin 2^0 + ^ sm 2.po + 2ro^ (^') J W^ 

ao2 = iro^'co" cos^ «9o + j- + "i (1 - 3 sin" ^0). 
"0 "0 
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On substituting (12) in the first of (13) and equating coefficients of corre- 
sponding powers of h, it is found that 

,^ ., — rpdi'-''^ [aio^Oo2 — auOoiaio + Ooi^azo] 

^^*^ "^~(a2<''))2 + ro^(62W)»' "' "= aio« 

On substituting (12) in the second and third equations of (6), it is found 
that 

^^ = ^» ^ (^T^OM^^TiW^'' "^ L"^\ dA" jo + "^^^ 
(16) +aiV»>]/i^+ •••, 

5. Determination of the Coordinates of P2. The point P2 is determined 
by the conditions that it shall he on the surface of the rotating body and that 
the normal to the equipotential surface passing through it shall also pass 
through Po' The first of these conditions gives an equation which differs 
from (8) only in that the subscript 2 appears in place of 1. 

Let 
(16) W = |wV cos^i <p + V. 

Then the equations of the normal to the eqvupotential surface TF = C at 
the point P2(a;2, 2/2, 22) are 



m) (m m) 

\ dx )i \dy Ji \ dz )i 



(17) \ bx Ji _ \ dy /2 _ \ dz J2 

X — X2 y ~ Vi z — Zi ' 
If this fine passes through Po these equations are satisfied by 
X = Xo = ro(l + h) cos (po cos do, 

(18) y = yo = roil + h) cos ^0 sin ^0, 
z = Zo = ro(l + h) sin (f>o. 

The two equations (17) and the one which expresses the condition that P2 
shall be on the surface of the rotating body are sufficient to determine 
Xi, 2/2, and Zi, or preferably their equivalents in polar coordinates. 

Suppose 7 is a sum of zonal harmonics. Then the two equations (17) 
can be replaced by a single one, which is, if ^0 is taken equal to zero, 
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The equation which expresses the condition that Pj shall be on the surface 
W = C and this equation become respectively in polar coordinates, after 
some simplifications, 

Fi = Wi^i^ cos'' <pi — To^ cos" ^o) + a ( - — - j 



(19) 



G^^ 



r-od + h) 



+ K^-^)-^K 



sm^ ^2 sm^ ^0 



n" 



W 



) = o, 



I « + ^ (1 - 5 sin" ^2) J sin (^2 - <(>o) 



6^, 



H — 4 [»'o(l + h) cos <po — Ti cos v'2] sin ^52 

+ r2w"[ro(l + ^) sin <p(, — r^ sin (p^ cos ^2 = 0, 



from which r^ and ^52 can be determined. 

In order to solve (19) it is convenient to let 



(20) 



r2 = ro(l + p), <Pi = <Po-\- a. 



Then it follows that F2 and (?2 can be expanded as power series in p, a, and /i 

of the form 

(21) F2 = /2(p, a, h) = 0, G2 = P2(p, <r, /i) = 0. 

It is noted that /2 does not involve h, g^ contains h Unearly, and both f^ 
and gi converge for any finite values of a and h provided \ p\ < 1. 

Both of equations (21) are satisfied by p = a = h = 0. The con- 
dition that they shall be imiquely solvable for p and a as power series in h, 
vanishing for h = 0,is that the Jacobian of /2 and g^ with respect to p and 
<r shall be distinct from zero for p = a = h = 0. It is found from (1ft) that 



(22) 





d/2 a/2 




/ = 


dp' da 

dgi dgi 
dp' da 


= 



2roa2(«, +2ro"b2<''^ 
2ro"&2<•^ -2ro02»> 



= - 4ro"[(a2«»)" + ro"(62<''>)"]. 



which is distinct from zero because it is — r^ times the square of the 
acceleration of gravity at Po<"^ Therefore the solution of (21) has the 
form 

(23) P = )3i/i + ^2/i"+ •••, (r = 71/1 + 72/1' + •••, 

which converge provided | fe | is sufficiently small. 

In order to determine the coefficients of (23) it is necessary to have the 
expUcit expansions of (21). On referring to (19) it is seen that these series 
have the form 
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Fi = 6100P + ftoioff + &200P^ + &iioP<r + fcojoo"^ + • • • = 0, 

(24) G2 = Cioop + CoioO" + Cooife + C200P* + Ciiopo- 

+ CioipA + Co2o<'^ + Coiio-A + . • . = 0, 

where all terms which are not zero up to the second order in p, a, and h 
have been written. It is easily found from equations (19) that the explicit 
values of the coefficients of (24) are 

bioo = 2ro02<"», &010 = 2ro2&2<'*>, 

On 

6200 = ko^'w" cos^ ^0 + r + :r3 (1 ~ 3 sm'' ^0), 

'0 ~0 

&110 = I - '"o^'w" + ^Ti I sin 2^90, 

&020 = - ko^w" + ^ cos 2^0, 
Cioo = 2u^hi'-'>\ Coio = - 2roa2(»>, 
^^^^ Cooi = - 2ro^62(«, C200 = [- W<^'' + ^] sin 2*^0, 
cuo ro''coH2 - 3 smVo) - ^ - ^ (3 - 8 sin* ^0), 

'0 '0 

Cioi = I |j-o*w'' - -^ I sin 2^50, 

C020 = I fro^w* - Y ;^3 J sin 2^o, 

Con = - r,?u>^ sinVo + J- + ri (3 - 7 sin* ^0). 

"0 '0 

In order to obtain the meridional deviation it is sufficient to compute 
71 and 72, but in order to obtain 72 it is necessary to compute |8i. It is 
foimd on substituting equations (23) in (24) and equating to zero the 
coefficients of corresponding powers of h, and then solving the resulting 
linear equations, that 

. ro*(b2<")' 

'^^ ~ (a2»))* + roH&2»0" 

(26) 'I'l - (a2("»)* + ro*(&2»')" 

■D72 = — &iooCoii7i — (&100C200 — Cioo&2oo)|8i'' 

— (&100C110 — Cioo&iio)^i7i ~ (&100C020 ~ Cioo&02o)7l* ~ &100Cl01^1, 
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where 

(27) D= - 4ron(c2«»>)^ + uKhi'-^m- 

These results substituted in (23) and (20) give the coordinates of Pi when 
F is a sum of zonal harmonics. 

6. Expressions for the Deviations. In case F is a sum of zonal har- 
monics Xj is zero, and it follows from the second equation of (15) that the 
expression for the angular deviation in longitude is 



(28) Xi = { 



— raOi 



(0) 



(02«>)^ + ro^(62«»)' 



C3»>/l* + 



The value of h, which is entirely arbitrary, can be taken so small that the 
sign of X will be determined by the first term of this series. 
It follows from (7) and (2) that 

(29 ^'' ^'^ 

62<"^ = — iw^sm 2^0 — 2~5si'i ^fa- 

The quantity 02'*^ is always negative for an actual physical body for other- 
wise there would be an acceleration outward along the radius at Po*"'- In 
the case of the earth the first term of the expression for Oj*"' is mmaerically 
about 300 times the coefficient of cos^'^o in the second term, and about 
600 times the coefficient of (1 — 3 sin^ ^0) in the third term. That is, 

(30) -» = SOOrow^ = 600 —: 

approximately. Consequently 

C3(« = -|«[^-b2Wtan<.o] 
(31) 

is positive. Therefore the deviation is toward the eastward. Since 62^"* 
is small compared to 02<''\ an approximate expression for the easterly 
deviation, obtained by neglecting bj'"' in the denominator of (28) and the 
terms in 02^°' and Cz'-"^ which are multiplied by w^ or )3, is 

(32) Xi = (^^fa,+ -.-);i?+ •••. 
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By making use of (23) and the explicit value of 71 given in (26), the 
meridional deviation is found from (15) and (20) to be 

(33) <pi - <P2 = Oh + ^ai(~^^ + a262»> + «iV« -72^'+ •••. 

The first important fact to be noted is that the coefficient of fe is identically- 
zero. It follows from the first equation of (15) and from (29) and (30) 
that the numerator and denominator of the coefficient of /i in the expression 
for <pi are polynomials in w* and /3, and that the quotient can be expanded 
as an infinite converging power series in w^ and /3. Exactly the same is of 
course true for the coefficient of /i in the expression for (pi. The quantities 
(Pi and <P2 are derived by quite different processes, the former involAong the 
solution of differential equations and the latter only the solution of implicit 
fimctions. Hence it is clear that methods of approximation in canying 
out these different processes are beset with danger because the effects upon 
the results are apt not to be the same id both. 

The coefficient of k^ La (33) is exact, but ia order to determine its sign 
all terms except those of lowest order in w^ and /3 may be neglected. It 
follows from the numerical coefficients which are involved and the relations 
given in (30) that, at least for the earth, these simpMcations can not 
change the sign of the result. It follows from (29) that 02^"' is of order zero 
in w^ and |8, and that bj^**' is of order one in these same quantities. Then 
it is seen from (7) that 64^°' and {db2i'dh)o are of order one in w" and p, though 
most of the terms in bi'-"^ are of higher order. They are explicitly 



(34) 



64(0) = _ — — ^ I 4^,2 -I j sin 2<po + terms of the second order, 

/db2\ ^150 . - 



It follows from (14) that ai and a2 are of order zero in w^ and 0, and from 
the first two equations of (26) that /3i is of the second order and 71 is of the 
first order. 

Now consider 72. Since D is of order zero all the terms in the right 
member of the third equation of (26), except possibly the first, are of the 
second order at least. It is seen from the first and last equations of (25) 
that bioo and Cou are of order zero. Therefore the first term of 72 is of order 
one. 

On retaining only the terms of the first order in the coefficient of ¥ in 
(33), it is found that the approximate expression for the meridional devi- 
ation is 
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(35) <.i-<.2 = [|-^(^-4«^)sin2<.o+ •••]/i'+ •••• 

For the eaxth Au? is more than five tunes as great as 5/3/ro^. Therefore the 

deviation of a body freely falling a small distance near the earth's surface is 

equatorward for all latitudes between and =*= 90°. 

The University op Chicago, 
December 1, 1913. 



